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INTRODUCTION 



Larainarization of the boundary layer on the surface of aircraft wings can be 
accomplished by the use of concepts such as Natural Laminar Flow (NLF) , Laminar-Flow 
Control (LFC) which uses suction over the entire surface of the wing, and Hybrid 
Laminar-Flow Control (HLFC) which uses suction only over 15 to 20 percent of the wing 
upper surface near the leading edge. The experimental data obtained at NASA Langley 
for the LFC wing have demonstrated that substantial reductions in wing profile drag 
can be obtained at transonic Mach numbers by the use of suction over the entire sur- 
face of the airfoil. This drag reduction is primarily the result of the presence of 
laminar boundary-layer flow over a large portion of the airfoil surface. . 

Extension regions of laminar flow can also be maintained on wing surfaces with 
the NLF concept, which involves the appropriate choice of pressure distribution to I 

limit the amplification of the disturbances that trigger the transition of the lam- 
inar boundary layer to a turbulent boundary layer. Flight tests at the Dryden Flight 
Research Center on a variable sweep TACT F — 111 fighter aircraft with an NLF wing 
glove have shown that laminar flow can be maintained over a large portion of the wing 
surface at transonic speeds. However, NLF is restricted to flight applications at j 

low Reynolds number conditions and to wings with relatively low sweep angles. j 

HLFC which combines the features of both LFC and NLF applies suction forward of j 

the front wing— box spar to prevent the transition of the laminar boundary layer due 
to cross— flow and attachment - line instabilities that occur on swept wings. Laminar j 

boundary-layer flow is maintained aft of the front spar by the selection of a pres- j 

sure distribution which surpresses the growth of the disturbances due to Tollmien- i 

Schlichting waves or due to cross flow. Although the extent of laminar flow is less 5 

for HLFC than for LFC, the conventional wing box structure can be retained and the j 

mechanical complexities are not as great. 

The design of the HLFC wing sections to minimize the wake drag requires the ; 

optimization of the drag due to laminar skin-friction and the turbulent boundary- 
layer separation. The trailing— edge recovery pressure is fixed due to the Kutta 
condition for subsonic and transonic flow and the aft pressure gradient, free-stream I 

Reynolds number, and forward suction levels determine whether or not the turbulent - 

flow will remain attached close to the trailing edge. Examination of existing exper - 

i mental data on the LFC wing has shown that if turbulent separation occurs upstream f 

of the 95-percent chord position on either the upper or lower wing surfaces, a pres- f 

sure distribution necessary to maintain laminar flow cannot be realized. Therefore, g 

the design of HLFC wings sections necessitates the use of reliable theoretical ; 

methods which accurately predict the locations of laminar boundary-layer transition I 

and turbulent boundary-layer separation. S 

Several finite— difference boundary layer, stability, and full Navier— Stokes 
c qua t ion solvers are available and have produced very encouraging results. However, I 
these methods are not well suited for the routine optimization studies that were 
performed during the design of the HLFC wing section. Several new integral | 

boundary-layer methods, which are applicable to swept wings with varying amounts 
of surface suction, have been developed for the prediction of laminar, transition, 
and separating turbulent boundary layers. These methods have been developed for use 
at either subsonic or supersonic speeds, have small computer execution times, and 
are simple to use. The purpose of this presentation is to briefly outline the 
theoretical equations and assumptions which form the basis of these boundary-layer 
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methods and to present the results of several correlation cases with existing experi 
mental data. The results of the application of these methods to the design of the 
HLFC wing scheduled to be tested in the Langley 8-Foot Transonic Pressure Tunnel will 
also be presented. 
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THEORETICAL METHODS USED DURING DESIGN OF HLFC WING 


The theoretical methods used during the design and analysis of the HLFC wing are 
listed in figure 1 and are classified as 1) potential/viscous design and analysis, 2) 
boundary-layer analysis, 3) wake and suction drag computations, and 4) a modified 
strip method for finite wings. 

Design of the upper surface of the HLFC wing in the local supersonic region at 
the design transonic free-stream Mach number was performed using the Perturbation 
Method of Characteristics technique* Modifications to the HLFC airfoi^ section were 
analyzed at transonic conditions using the Korn/Garabedian airfoil code and at sub- 
sonic conditions using the Multi-Component (MCARF) airfoil code. 

Analysis of design modifications on the viscous characteristics of the HLFC wing 
were performed using several recently developed integral boundary-layer methods* 

These methods consist of 1) an Integral compressible laminar boundary-layer method 
for swept wings in the presence of suction at subsonic and supersonic speeds, 

2) criteria for prediction of laminar boundary-layer separation and reattachment, 

3) criteria for prediction of the location transition due to either the amplification 
of Tollmien-Schlichting waves, cross flow, oi^leading-ed ge contamination, and 4) new 
integral separating turbulent boundary layer* 

Theoretical methods were also developed to determine the effect of design modi- 
fications on the drag characteristics of the HLFC wing. These methods account for 
both the changes in the wake drag and suction drag as a result of applying suction on 
the upper surface leading-edge region of the HLFC airfoil. 

A modified strip method was developed during this design study to account for 
the finite and swept wing properties of the HLFC wing* This strip method also 
accounts for taper and both spanwise and chordwise pressure gradients. 


• Potential/viscous design and analysis 

• Perturbation method of characteristics for inverse design at transonic 
speeds 

• BGK IBauer-Garabedian-Korn) for transonic analysis 

• MCARF for subsonic analysis 

• Boundary-layer analysis 

• Integral compressible laminar boundary layer with sweep and suction at 
subsonic thru supersonic speeds 

• Short bubble and reattachment criteria 

• T.S. and C.F. transition criteria 

• Separating turbulent boundary layer method (AIAA 86-1832-CP) 

• C n „ and C n 

U, wake U, suction 

• Modified strip method 


Figure 1 


550 



GOVERNING EQUATIONS 


The usual governing equations for compressible hydrodynamic laminar boundary- 
layer flow which consist of the continuity, streamwise-momentum, normal-momentum, and 
cross-flow-momentum equations are presented in figure 2. These equations contain the 
terms consisting of variable physical properties, such as, density, p, and dynamic 
viscosity, y. The values of these physical properties vary across the boundary layer 
as well as along the flow direction, and these variations are non-negligible for 
boundary-layer flow at transonic speeds* At supersonic speeds, these variations in 
the physical properties of fluid within the boundary layer are quite appreciable. 

This means that there is a strong coupling between solutions of hydrodynamic and 
thermal boundary-layer equations at high transonic and supersonic speeds. 

In order to simplify the governing equations for solution by Integral techniques 
while maintaining realistic, computational results for the hydrodynamic and thermal 
boundary layers, Stewartson's transformations are used. These transformations refor- 
mulate the boundary-layer equations of motions into a transformed plane which is 
Independent of the varying physical properties of fluid. The relations between the 
velocities in the transformed and physical planes are also shown In figure 2. 
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TRANSFORMED BOUNDARY-LAYER EQUATIONS AND BOUNDARY CONDITIONS 

The transformed boundary-layer equations using Stewartson's tranf ormations are 
presented in figure 3. These tranforraed boundary-layer equations are applicable to 
infinite swept wings. The relation between the pressure gradient in the transformed 
and physical planes and the transformed boundary conditions needed to derive the in- 
tegral equations is also shown. The several groups of physical dimensionless param- 
eters used in the development of the integral compressible boundary-layer method are 
also presented. The subscript n w" indicates that the parameter is defined at the 
wall. 
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COMPATIBILITY CONDITIONS AND STREAMWISE MOMENTUM-INTEGRAL EQUATION 

WITH SUCTION 


The compatibility conditions used to derive the integral boundary-layer equa- 
tions are presented in figure 4. The streamwise momentum-integral equation with 
suction included was derived by integrating the transformed boundary-layer equations 
from the wall to the edge of the boundary layer and by making use of Leibnitz’s rule. 
Furthermore, by making use of the dimensionless parameters defined on figure 3, a set 
of simultaneous equations is derived for the solution for the boundary-layer momentum 
thickness and form factor in the transformed plane in the presence of suction. The 
form of function F(S,K g ) is derived by the curve fit of the exact solution results 
for the Falkner-Skan type flow in the presence of suction. 
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Figure 4 
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BOUNDARY-LAYER PARAMETERS AND CROSS FLOW MOMENTUM-INTEGRAL EQUATION 


The application of reverse Stewartson's transformations yields the relationships 
between the parameters in the physical and transformed planes as presented in figure 
5. The transformed cross-flow momentum-integral equation is derived in a manner sim- 
ilar to that for the streamwise momentum-integral equation. The numerical solution 
of this equation gives the computational results in the transformed plane in the form 
of the ratio of the boundary-layer thickness in the cross-flow direction to that in 
the streamwise direction as a function of the dimensionless parameter, S, which is 
the ratio of the distance along the normal section in the streamwise direction to the 
chord of normal section. The cross-flow boundary-layer thickness, 6 , can be calcu- 
lated by multiplying the ratio 6 z /6 g , which is obtained by solving the cross-flow 
momentum-integral equation, by the streamwise boundary-layer thickness 6 g . The 
dimensionless shape of the velocity profile w/W g can immediately be calculated from 
the knowledge of <5 Z , momentum thickness, and displacement thickness. 
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THERMAL BOUNDARY-LAYER EQUATIONS 


Equations for thermal boundary layers in transformed and physical planes are 
presented in figure 6. As a first approximation, assuming that the external pressure 
gradient is negligible (which is approximately the case for the supersonic wing at 
low lift cruise conditions) and that the Frandtl number is close to unity, an analyt- 
ical expression for the temperature profile can be derived* This expression is then 
modified as shown to account for real pressure gradients and the exact value of 
Prandtl number. The real pressure gradient is accounted for through the polynominal 
expression for the velocity ratio u/U g where the A^ T s are functions of the local 
pressure gradient* Once the temperature profiles are computed, the variations of 
dynamic viscosity, kinematic viscosity, and density across the boundary layer can be 
calculated and their average values determined* These averaged values are then 
utilized to compute dimensionless parameters which are used to determine the transi- 
tion location of the laminar boundary layer. 
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CRITERIA FOR DETERMINING LOCATION OF BOUNDARY-LAYER NEUTRAL INSTABILITY 


The relationship between the dimensionless pressure-gradient parameter and an 
"equivalent 11 Reynolds number based on momentum thickness for neutral instability is 
presented in figure 7* The curve presented Is used for determining the location of 
neutral instability for either the streamwise or cross-flow laminar boundary layer. 
The effects of pressure gradients, Mach number less than 1.3, and suction were 
accounted for during the derivation of this curve from the solution of Orr-Sommerfeld 
equations in conjunction with Stewartson’s transformations. 
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EFFECT OF FREE-STREAM TURBULENCE ON DIFFERENCE BETWEEN TRANSITION- 
AND NEUTRAL-INSTABILITY REYNOLDS NUMBER FOR FLOW ON A FLAT PLATE 

The effect of free-stream turbulence on the transition location on a flat-plate 
laminar boundary-layer subjected to zero pressure gradient is shown in figure 8. The 
transition data were obtained by several investigators in the range of free-stream 
Mach numbers from low subsonic to a supersonic value of 3. The temperature profile 
data used to construct this curve were obtained either from experimental data or cal- 
culated from the previously stated theoretical expression for the temperature profile 
presented in figure 6. The curve shows the difference in "equivalent” Reynolds num- 
ber at the transition location and at the point of neutral stability. The equations 
for the averaged value of kinematic viscosity v and "equivalent" integral thickness 
6* which were used to derive this curve from the experimental data are also 
presented. 


^tran ^inst 


U gbn 9 I I T 

and t, -0. ) T d 1 



Figure 8 
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BOUNDARY-LAYER TRANSITION CRITERIA 

The criteria used in the compressible laminar boundary-layer method to determine 
the transition location due either to the amplification of Tollmien— Schlichting waves 
or to cross flow is presented in figure 9. The dimensionless expression on the ab- 
scissa contains several implicit and explicit physical parameters that are signifi- 
cant during the transition of laminar boundary layer. In deriving the curves shown 
in this figure for several values of free-stream turbulence intensities, use has been 
made of information presented in figure 8. It was assumed that the effect of free- 
stream turbulence on the transition of laminar boundary layer without a pressure gra 
dient behaves in a manner similar to a laminar boundary layer with a pressure 
gradient. 



Figure 9 
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EQUATIONS FROM APPLICATION OF SWEEP THEORY 


A schematic illustrating the application of "Strip Theory" for the prediction of 
laminar boundary-layer transition on a finite swept, tapered, and twisted wing is 
presented in figure 10. The wing is divided into a large number of finite strips 
oriented normal to the leading edge (more accurately normal to the local aerodynamic 
center-of-pressure line)* The pressure distribution Cp on the edge as well as on the 
centerline of each strip can be determined either by the available two- or three- 
dimensional theoretical flow codes. The local Mach number distributions In the 
direction parallel and normal to the strip are also presented in this figure. These 
local Mach number distributions, Mj and n , are Input Into the compressible lam- 
inar boundary-layer method to determine whether or not transition will occur due to 
either Tollmien-Schlichting waves or cross-flow instabilities. 
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Figure 10 
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CORRELATIVE STUDIES 


The present Integral boundary-layer method for predicting transition was used 
extensively during the design of the Hybrid Laminar-Flow Control (HLFC) wing to be 
tested in the Langley 8-Foot Transonic Pressure Wind Tunnel. In order to establish 
the validity of the method prior to design of the HLFC wing, the correlative studies 
listed in figure 11 were performed to compare the computational results by the pre- 
sent theory with available experimental data. 


• Variable sweep wing with NACA 642A015 section 


• Phoenix wing 


• NASA LFC wing 


Figure 11 



EXPERIMENTAL SETUP FOR TRANSITION MEASUREMENTS ON NASA AMES VARIABLE SWEEP WING 


8 

A schematic of the experimental setup used by Boltz in the NASA Ames Research 
Center’s 12-Foot Pressure Tunnel to determine the transition location on an infi- 
nitely swept wing at various sweep angles is shown in figure 12* The sweep angle was 
varied from 0° to 50° by insertion of wedges, and the wing tips were kept parallel to 
the free stream with appropriate wing tip extension. The wing section was an NACA 
642A015 airfoil. The chord of the wing was four feet, and the wing was mounted ver- 
tically on the turntable in a semispan manner. In the unswept position the wing had 
an effective aspect ratio of 5.0. The transition locations were experimentally de- 
termined at the various sweep angles and Reynolds numbers using both flow visualiza- 
tion techniques and small microphones located in the model. 





COMPARISON OF PRESSURE DISTRIBUTION FOR NACA 64 2 A015 WING WITH 30° SWEEP 

The pressure distribution computed by the NASA Multi— Component Airfoil Program 
(MCARF) for the NACA 64nA015 section at an angle of attack of 1° and sweep angle of 
30° is shown in figure 13. The experimental data shown in this figure were obtained 
by Boltz at the same free-stream conditions and sweep* The comparison shows excel- 
lent agreement between theory and experiment* 


A = 30°, a = 1°, M co = 0.2 7 



Figure 13 
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COMPARISON OF PRESSURE DISTRIBUTION FOR NACA 64 2 A015 WING WITH 50° SWEEP 

The comparison between the theoretically computed pressure distribution and 
experimental data by Boltz for an angle of attack of 0° and a sweep angle of 50° is 
shown in figure 14. The computed pressure distribution agrees quite well with the 
experimental data except In the vicinity of the trailing edge* This discrepancy at 
the trailing edge can be attributed to turbulent separation due to outboard washout 
phenomenon which usually occurs at this relatively high sweep angle* 


A = 50°, a=0°, Moo = 0.27 



Figure 14 
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TRANSITION CORRELATION FOR NACA 64 2 A015 WING WITH 20° SWEEP 

The comparison between predicted and experimentally measured transition loca 
tions for the NACA 64 2 A015 wing with 20° of sweep Is presented in figure 15. The 
theoretical results are shown as a function of free-stream Reynolds number for tran- 
sition locations due to 1) Tollmien-Schlichting waves or laminar separation with 
short bubble turbulent reattachment and 2) transition due to cross flow. Also shown 
in this figure are the computed x/c locations for Tollmien-Schlichting (T.S.) and 
cross-flow (C.F.) neutral Instabilities as a function of free-stream Reynolds number. 
The correlation results suggest the following: 

(1) Predicted transition due to T.S. by the present method agrees quite well 

with experimental data for Reynolds numbers less than 20 million. 

(2) For Reynolds numbers greater than 20 million, the predicted transition due 

to C.F. occurs upstream of that due to T.S., and the predicted locations 

due to C.F. agree well with the experimental data. 

(3) For a sweep of 20° and Reynolds numbers greater than 20 million, the 

transition is due completely to C.F. 
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Figure 15 



TRANSITION CORRELATION FOR NACA 64 2 A015 WITH THE 40° SWEEP 

The correlation of results between present method and Boltz's experimental data 
for the higher sweep angle of 40° is presented in figure 16. The theoretical data 
presented were computed using the pressure distribution generated by the MCARF pro- 
gram. The following conclusions can be drawn from the correlations presented in this 
figure: 

(1) Transition locations due to T.S. and C.F. are predicted fairly well. 

(2) For Reynolds numbers greater than 8 million, theoretical computations indi- 

cate that C.F. triggers the transition more abruptly and earlier than at 
the lower sweep angles. 

(3) There is a danger of leading-edge contamination due to C.F. for Reynolds 

numbers larger than 18 million at a sweep angle of 40°. 


A = 40° ( a = 0°, Mjo'0.27 


x/c 

transition 



Figure 16 
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PLANFORM OF PHOENIX SAILPLANE AND AIRFOIL GEOMETRY 

The planform of the Phoenix sailplane and the geometry of the airfoil at the 
test location on the wing are shown In figure 17, The pressure distributions and 
velocity profiles were measured by Raspet^ at several chordwise locations at the test 
location for several values of lift coefficient* The maximum thickness-to-chord ra- 
tio for the test location airfoil Is 0.15 and is located at x/c *= 0.35. The maximum 
camber is located at x/c * 0.70. The aspect ratio of the wing is approximately 18. 


Max thickness 14% Max camber 4.5% 





PRESSURE DISTRIBUTION CORRELATION ON PHOENIX WING 
C L = 0.76 AND R n = 1.464 x 10 6 

The comparison between the theoretical pressure distribution computed by the 
MCARF program and the experimental data of Raspet for the Phoenix wing at a lift 
coefficient of 0.76 and a chord Reynolds number of 1.464 million is shown in figure 
18. The results plotted in this figure show that the computed pressure distribution 
3grees quite well with the measured experimental data. These computed and experl 
mental pressure distributions were Input as boundary conditions to the integral 
boundary-layer method used to compute the transition locations presented 
subsequently. 


Cl=0.76, R n = 1.464 X 10 6 



Figure 18 
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PRESSURE DISTRIBUTION CORRELATION ON PHOENIX WING 
C L = 1.14 AND R N = 1.2 x 10 6 

The comparison of pressure distribution computed by MCARF and the experimental 
data for Phoenix wing at lift coefficient of 1*14 and a chord Reynolds number of 
1.2 million is presented in figure 19. The computed pressure distribution agrees 
fairly well with experimental data except near the trailing edge where turbulent 
boundary-layer separation is present. 


i 


C L = 1.14, Rn * 1.2 x io 6 



Figure 19 
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TRANSITION CORRELATION FOR PHOENIX SAILPLANE WING 


The correlation for transition and turbulent separation locations on the test 
location airfoil of the Phoenix wing is presented in figure 20. The results are 
plotted as a function of lift coefficient. These results indicate that the chordwise 
location of turbulent separation is predicted fairly well. However, some discrepancy 
exists between the computed and experimental transition location for lift coeffi- 
cients in the range of 0.6 to 1.1 when the theoretical pressure distributions were 
used as input boundary conditions. In order to examine the effect of using different 
boundary conditions, transition locations were also computed using experimental pres- 
sure distributions as input boundary conditions. As seen In figure 20, the differ- 
ence in the theoretical transition locations using the two different boundary condi- 
tions Is not significant. It should be emphasized that the author Raspet did not 
measure the location of transition directly but, instead, inferred the location from 
the measured velocity profiles; therefore, there may be some discrepancy in the 
actual measured transition locations. 
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Turbulent 
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(1) Computed transition with B.C. of computed pressure distribution 

(2) Computed transition with B.C. of experiment pressure distribution 

(3) Computed turbulent separation with B.C. of computed 



Figure 20 
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CORRELATION OF TRANSITION DATA FOR LFC WING WITH PARTIAL SUCTION 

A = 23° AND = 0.826 

The correlation of the transition locations for the NASA Langley LFC wing with 
partial suction on the upper surface are presented in figure 21. The chordwise 
location of suction was varied from an x/c of 0.0 to 0.5. The transition measure- 
ments were made at a free-stream Mach number of 0.826 and chord Reynolds numbers of 
10 to 20 million. Theoretical computations were performed using boundary conditions 
of experimentally measured suction coefficients as input boundary conditions. The 
predicted transition locations are in good agreement with the measured transition 
locations. The results of the three correlation cases presented have demonstrated 
the validity of the new integral boundary-layer method. 


x/c 

transition 



Figure 21 



DESIGN CONSTRAINTS OF HYBRID LAMINAR-FLOW CONTROL WING 


The design constraints and objectives for the Hybrid Laminar Flow Control (HLFC) 
wing to be tested in the Langley 8-Foot Transonic Pressure Tunnel are illustrated on 
figure 22. The design constraints and objectives that were laid out by Mr. P. J. 
Bobbitt, who also directed the HLFC design and optimization studies, are summarized 
as follows: 


Design Constraints 


(1) The geometry of the lower surface of the HLFC wing must be the same as that of 

the LFC wing. 

(2) The geometry of the upper surface panel #1 of the HLFC wing must be identical 

to that on the present LFC wing and suction must be applied only through 
panel #1. 

(3) The sonic bubble height for the HLFC wing must be no greater, and, if possi- 

ble, smaller than that for the LFC wing. 


Design Objectives 


(1) The geometric shape of the upper surface panels #2 and #3 must be derived by 

the inverse perturbation method of characteristics so that laminar boundary- 
layer flow is maintained up to x/c = 0.6 on the upper surface of HLFC wing 
for C L = 0.45 and M^ = 0.82. 

(2) The computed turbulent separation location on the upper surface of the HLFC 

wing must be aft of x/c = 0.95. 



• Suction applied to panel 1 

• Geometry of panels 2 and 3 altered for: 

• Sonic line height constraint 

• Laminar B.L. flow for x/c % 0.6 

• Turbulent separation aft of x/c = 0.95 

Figure 22 
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COMPARISON BETWEEN THEORETICAL AND EXPERIMENTAL PRESSURE DISTRIBUTIONS 

FOR LFC AND HLFC WINGS 


A comparison of pressure distributions between the design HLFC wing and the pre- 
sent LFC wing is presented in figure 23. The application of suction is limited to 
0.025 < x/c < 0.26 on the upper surface of the HLFC wing. The pressure distribution 
for the HLFC wing (shown by the dotted line) constitutes the boundary condition 
necessary to achieve the design constraints and objectives outlined on figure 22. 



Figure 23 
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PREDICTED TRANSITION AND TURBULENT SEPARATION LOCATIONS FOR HLFC 


The computational results for the predicted transition and turbulent boundary- 
layer separation locations on the upper surface of the HLFC wing are presented in 
figure 24. Predicted results are shown for transition due to Tollmien-Schlichting 
(T.S.) wave amplification, transition due to cross flow (C.F.), and instability due 
to T.S. and C.F. The chordwise location of computed turbulent separation is also 
presented. These results are plotted as a function of free-stream Reynolds number 
and for a wing C^ * 0.44 and free-stream = 0.811. 

As shown in figure 24, transition due to T.S. with short bubble reattachment 
occurs at x/c = 0.58 for the range of Reynolds numbers shown. However, transition 
due to C.F. takes over at a Reynolds number of 15 million. It is assumed with the 
present theory that the transition occurs on the wing due to whichever phenomena 
appears first. Thus, cross flow essentially determines the transition location at 
the higher Reynolds number. In addition, there is a danger of leading-edge contami- 
nation due to C.F. instability at Reynolds numbers larger than approximately 20 mil- 
lion. The computed turbulent separation location is downstream of x/c = 0.95 for the 
entire range of Reynolds numbers. 


A = 23°, Cl = 0.438, Moo = 0.811 



Figure 24 
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CONCLUSIONS AND RECOMMENDATIONS 


The results of the design studies of the HLFC wing using the integral boundary-layer 
methods of Goradia lead to the following conclusions and recommendations: 

(1) The perturbation method of characteristics was found to be extremely useful for 
determining the geometric changes needed in the baseline LFC airfoil shape to obtain 
the desired HLFC pressure distribution. 

(2) The compressible laminar boundary-layer method with suction was found to be 
quite accurate in predicting the extent of laminar flow for swept wings in the pre 
sence of suction. Also with this method the suction requirements can be determined 
for different local velocity profiles such as those for Blasius, separating laminar, 
or asymptotic suction velocity distributions. 

I 

(3) The turbulent boundary-layer separation method includes terms to account for the ! 

rapid increase in the turbulent fluctuations in the flow near separation; these terms 
greatly improved the accuracy of the prediction of the location of turbulent : 

separation. 

(4) The integral boundary-layer methods execute very rapidly on the computer making | 

it possible to analyze several hundred configurations in a relatively short period of 

t ime . | 

(5) The results of the correlative studies generally showed excellent agreement be- 
tween the theoretical predictions and the experimental data. The results also showed j 

that, for values of Reynolds number and wing sweep of practical interest for commer— ; 

cial and fighter aircraft, cross-flow instabilities were predominant in triggering ! 

transition. ! 

(6) The analysis of the final HLFC design showed that less than one count of suction 
drag coefficient was required in conjunction with the appropriate pressure distribu- 
tion to achieve laminar flow on the upper surface of the wing to the 60-percent chord 
location. This one count of suction drag results in a corresponding 30 to 40 count 
reduction in the wake drag coefficient which suggests that HLFC is a very lucrative 

and promising concept for viscous drag reduction at both transonic and supersonic . 

speeds. . 

(7) In order to prevent laminar transition due to cross flow at large values of 

sweep and Reynolds number, it is recommended that both the chordwise and the spanwise 
pressure gradients be tailored to minimize the growth of the boundary-layer distur- I 

bances. With the use of a modified strip theory, arbitrary spanwise pressure gradi- 
ents can be accounted for when using the present integral boundary-layer methods. 

The present integral boundary-layer methods have been extended to compute both lami- / 

nar and turbulent boundary-layer flow and to predict transition locations at super- : 

sonic Mach numbers. The effects of variations in temperature profile and physical 
flow properties across the boundary layer, suction, wing sweep, wing taper, and wing | 

are also accounted for in these methods . These methods are not limited to 
wings alone, but can also be used for the analysis of fighter aircraft fuselage with 
suction at supersonic flight conditions. Correlation studies are currently under way 
to determine the validity of these extended integral boundary-layer methods. jj 
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